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The shape of the hadronic form factor f+(q 2 ) in the decay D° — > K~e + v e has been measured 
in a model independent analysis and compared with theoretical calculations. We use 75 fb" 1 of 
data recorded by the BABAR detector at the PEPII electron-positron collider. The corresponding 
decay branching fraction, relative to the decay D — > K~tt + , has also been measured to be Rd = 
BR(D° -> K-e + v e )/BR(D° -> K~-k+) = 0.927 ± 0.007 ± 0.012. From these results, and using the 
present world average value for BR(D° — > K~tt + ), the normalization of the form factor at q 2 = 
is determined to be /+(0) = 0.727 ± 0.007 ± 0.005 ± 0.007 where the uncertainties are statistical, 
systematic, and from external inputs, respectively. 
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I. INTRODUCTION 

Measurements of exclusive semileptonic D decays pro- 
vide an accurate determination of the hadronic form fac- 
tors entering in these decays. Assuming that the CKM 
matrix is unitary, the elements \V CS \ and \ V c< i\ can be de- 
termined: 



\V a \ = \V ud \ - 



\V r 



cb 



0(A 6 ) = 0.9729 ±0.0003, (1) 



using the measured values [1] of \V u d\ and \V c b\, and the 
sine of the Cabibbo angle A = sin(# c ) ~ 0.227. Theo- 
retical predictions give estimates of the form factors in 
exclusive semileptonic B and D meson decays. Precise 
measurements of the hadronic form factors in D decays 
can help to validate predictions from QCD calculations in 
both D and B decays. Better understanding of the form 
factors in B decays is necessary to improve the precision 
on the determination of \V c b\ and \V u t\. 

In D° — > K~e + v e decays [2], with a pseudoscalar 
hadron emitted in the final state, and neglecting the elec- 
tron mass, the differential decay rate depends only on one 
form factor f + (q 2 ), 



dT 

dq 2 



Gf_ 

24tt 3 



\v cs \ 2 \p k (q 2 )\ \U(q 2 )[ 



(2) 



where Gf is the Fermi constant, q 2 is the invariant mass 
squared of the two leptons, e + and u e , and pxiq 2 ) is 
the kaon three- momentum in the D° rest frame [3]. In 
this paper we present measurements of the q 2 variation 
and absolute value of the hadronic form factor at q 2 = 
for the decay D° — ► K~e + is e (j). The data consist of 
D mesons produced in e + e~ — > cc continuum events at 
a center of mass energy near the T(45) mass, and were 
recorded by the BABAR detector at the Stanford Linear 
Accelerator Center's PEP-II collider. A semi-inclusive 
reconstruction technique is used to select charm semilep- 
tonic decays with high efficiency. As a result of this ap- 
proach, events with a photon radiated during the D° de- 
cay are included in the signal. The systematic uncertain- 
ties are kept as low as possible by using control samples 
extracted from data where possible. 

Measurements of D — > Klvi, based on smaller signal 
events samples, have been published by the CLEO [4], 
FOCUS [5] and Belle [6] collaborations.' 

This paper is organized as follows. A general descrip- 
tion of the hadronic form factor, f+(q 2 ), is given in Sec- 
tion II, where the different parameterizations considered 
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in this analysis are explained. In Section III a short de- 
scription of the detector components that are important 
to this measurement is given. The selection of signal 
events and the rejection of background are considered in 
Section IV. In Section V, the measured q 2 variation of 
the hadronic form factor is discussed and compared with 
previous measurements. In Section VI the measured de- 
cay rate is given and in Section VII these measurements 
are combined to obtain the value of /+ (0) . 



II. THE F+(Q 2 ) HADRONIC FORM FACTOR 

The amplitude for the decay D° — > K~t + vi depends 
on two hadronic form factors: 

<K(p')\V,\D(p)> = L +p ^-q^^l^k)f +{ q 2 ) 



+ 



m D- m K „ , / 2x 
~2 <W0W ) 



(3) 



where U M = s-f^c. The constraint /+(0) = /o(0) en- 
sures that there is no singularity at q 2 = 0. When the 
charged lepton is an electron, the contribution from / 
is proportional to ml and can be neglected in decay rate 
measurements. 

The parameterizations of f+(q 2 ) which have been com- 
pared with present measurements and a few examples of 
theoretical approaches, proposed to determine the val- 
ues of corresponding parameters, are considered in the 
following. 



A. Form factor parameterizations 

The most general expressions of the form factor f+(q 2 ) 
are analytic functions satisfying the dispersion relation: 



/+(<z 2 ) 



Res{f+) q 2 



1 



m 



dt- 



t — q 2 — i 



(4) 



The only singularities in the complex t = q 2 plane orig- 
inate from the interaction of the charm and the strange 
quarks in vector states. They are a pole, situated at the 
D* mass squared and a cut, along the positive real axis, 
starting at threshold (t+ = (mj) + mj() 2 ) for D°K~ pro- 
duction. 



1. Taylor expansion 

This cut i-plane can be mapped onto the open unit 
disk with center at t = to using the variable: 



z(t,t ) = 



to 



VU=t + - to 



(5) 



In this variable, the physical region for the semileptonic 
decay (0 < t < t_ = <^ ax = (mD — mx) 2 ) corresponds to 
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a real segment extending between ±z n 
value of z„ 



= ±0.051. This 
is obtained for to = t+ ^1 — \/l — t-/t+j . 
The z expansion of / + is thus expected to converge 
quickly. The most general parameterization [7], consis- 
tent with constraints from QCD, 



/+(*) 



1 



o fc (to) z k {t,t ), 



(6) 



is based on earlier considerations [8]. The function 
P(t) = z(t,m 2 D ,,) has a zero at the D* s pole mass and 
\P\ = 1 along the unit circle; $ is given by: 



*(Mo) 



24,„ V^)'(^ ri+ ^"" 



(^+-4+^+ - to) 

3 

+ - *-) 2 (*+ - ty 



(7) 



where xv can be obtained from dispersion relations using 
perturbative QCD and depends on u = m s /m c [9]. At 
leading order, with u — [10], 



Xv = 



-.2 ■ 



(8) 



with (5 polc = (l/7 po ie - a P oie)/(l - "polo) and /3 polc = 

l/7polo- 

If in addition, the form factors /+ and /o must obey a 
relation, valid at large recoil and in the heavy quark limit, 
then ctpoic = l/7 P oie [H] (/3 P oie = a P oic and 5 po \ c = in 
this case). Equation (10) becomes: 



Ma 2 ) 



/+(o) 



ttpolc 



(11) 



known as the modified pole ansatz. Initially an even 
simpler expression, the simple pole ansatz, was proposed 
which considered only the contribution from the D* pole. 
In the following, the pole mass entering in 



Mi 2 ) 



Mo) 



i - 



(12) 



is fitted. Note that such an effective pole mass value has 
no clear physical interpretation and that the proposed q 2 
variation does not comply with constraints from QCD. 
The obtained value may nonetheless be useful for com- 
parison with results from different experiments. 



The choice of P and $ is such that: 



5>*(to)<i- 



(9) 



k=0 



Having measured the first coefficients of this expansion, 
Eq. (9) can constrain the others. This constraint, which 
depends on xv, may have to be abandoned in the case of 
charm decays as the charm-quark mass may not be large 
enough to prevent the previous evaluation of xv from 
receiving large l/m c and QCD corrections. However the 
parameterization given in Eq. (6) remains valid and it 
has been compared [7] with available measurements. The 
first two terms in the expansion were sufficient to describe 
the data. 



2. Model- dependent parameterizations 



A less general approach assumes that the q 2 variation 
of f+{q ) is governed mainly by the D* pole and that 
the other contributions can be accounted for by adding 
another effective pole at a higher mass [11]: 



Mi 2 ) = 



MQ) 

1 - Ckpolc \ 1 



C^polc 



1 - 



= /+(0) 



3polo" 



1 - /3polo 



(10) 



B. Quantitative expectations 

Values of the parameters that determine f+(q 2 ) were 
obtained initially from constituent quark models and 
from QCD sum rules. These two approaches have an in- 
trinsically limited accuracy. In this respect, results from 
lattice QCD computations are more promising because 
their accuracy is mainly limited by available computing 
resources. 



1. Quark Models 

Quark model calculations estimate meson wave func- 
tions and use them to compute the matrix elements that 
appear in the hadronic current. There are a large variety 
of theoretical calculations [12]. Among these models we 
have selected the ISGW model [13], simply because it is 
widely used to simulate heavy hadron scmileptonic de- 
cays. This model was expected to be valid in the vicinity 
of Qmax) a region of maximum overlap between the ini- 
tial and final meson wave functions. In ISGW2 [14] the 
exponential q 2 dependence of the form factor has been 
replaced by another parameterization, with a dipole be- 
havior, expected to be valid over a larger q 2 range: 



/i SGW V) = 



/+(?max) 



(l + OI/(?max " ? 2 )) 2 ' 



ai 



12 



(13) 



The predicted values of the parameters are /+(<Zmax) = 
1.23 and r = 1.12 GeV" 1 [14]. 
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TABLE I: Parameterizations of f+(q 2 



modeling 



parameters 



expected values 



z expansion[8] do, r k — a k /ao 

general 2-poles [11] /+(0), /3 po ie, <5 P oio 

modified pole [11] /+(0), o? po ic 

simple pole /+(0), m po ie 

ISGW2[14] /+(*-), on 



no prediction 
no prediction 

<5 p ole = 
?™ P ole = m-Dj 

/+(«-) = 1.23 
ai = 0.104 GeV" 



2. QCD sum rules 

QCD sum rules [15] and their extension on the light 
cone [16], are expected to be valid at low q 2 . Using a 
value of 150 MeV for the strange quark mass, one obtains 
[16]: 

/+(0) = 0.78 ±0.11 and a polc - -0.071°;^, (14) 

using the modified pole ansatz. The uncertainty of /+(0) 
is estimated to be of order 15%, and the q 2 dependence 
is expected to be dominated by a single pole at the D* 
mass because the value of a po i c is compatible with zero. 



3. Lattice QCD 

Lattice QCD computation is the only approach able 
to compute f+{q 2 ) from first principles. Current results 
must be extrapolated to physical values of light quark 
masses and corrected for finite lattice size and discriti- 
zation effects. There have been several evaluations of 
f+(q 2 ) for different values of the momentum transfer in 
the quenched approximation [17, 18]. These results have 
been combined [17], giving /+(0) = 0.73 ± 0.07. The 
first unquenched calculation has been published recently 
[19]: /+(0) = 0.73 ± 0.03 ± 0.07 and a po i c = 0.50 ± 0.04, 
using the modified pole ansatz to parameterized the q 2 
dependence of the form factor. 



C. Analyzed parameterizations 

The different parameterizations of f+{q 2 ) considered in 
this analysis are summarized in Table I, along with their 
corresponding parameters and expected values, where 
available. 



with track elements in the 40-layer drift chamber (DCH) , 
which is filled with a gas mixture of helium and isobu- 
tane. Slow particles which do not leave enough hits in the 
DCH due to the bending in the 1.5-T magnetic field, are 
reconstructed in the SVT. Charged hadron identification 
is performed combining the measurements of the energy 
deposition in the SVT and in the DCH with the informa- 
tion from the Cherenkov detector (DIRC). Photons are 
detected and measured in the CsI(Tl) electro-magnetic 
calorimeter (EMC). Electrons are identified by the ra- 
tio of the track momentum to the associated energy de- 
posited in the EMC, the transverse profile of the shower, 
the energy loss in the DCH, and the Cherenkov angle in 
the DIRC. Muons are identified in the instrumented flux 
return, composed of resistive plate chambers interleaved 
with layers of steel and brass. 

The results presented here are obtained using a to- 
tal integrated luminosity of 75 ffcU 1 registered by the 
BABAR detector during the years 2000-2002. Monte Carlo 
(MC) simulation samples of T(4S) decays, charm and 
other light quark pairs from continuum equivalent, re- 
spectively, to 2.8, 1.2 and 0.7 times the data statistics, 
respectively, have been generated using GEANT4 [22]. 
These are used mainly to evaluate background compo- 
nents. Quark fragmentation, in continuum events, is de- 
scribed using the JETSET package [23]. The MC distri- 
butions have been rescaled to the data sample luminosity, 
using the expected cross sections of the different compo- 
nents (1.3nb force, 0.525 nb for B + B~ and B°B°, 2.09nb 
for light uu, dd and ss quark events). Dedicated sam- 
ples of pure signal events, equivalent to seven times the 
data statistics, are used to correct measurements for effi- 
ciency and finite resolution effects. They have been gen- 
erated using the modified pole parameterization ansatz 
for ,f+(q 2 ) with a po i c = 0.50. Radiative decays (D° — ► 
K-e+Vei) are modeled by PHOTOS [24]. To account for 
one of the most important sources of background, a spe- 
cial event sample with, in each event, at least one cascade 
decay D* + — > D°ir + , D° — ► K~Ti°e + v e (or its charge 
conjugate) has been generated with a parameterization 
of the form factors in agreement with measurements from 
the FOCUS collaboration [25]. Events with a D*+ and 
a D° decaying into K~tt + or K~ir + ir° have been recon- 
structed in data and simulation. These control samples 
have been used to adjust the c-quark fragmentation dis- 
tribution and the kinematic characteristics of particles 
accompanying the D meson in order to better match the 
data. They have been used also to measure the recon- 
struction accuracy on the missing neutrino momentum. 



III. THE BABAR DETECTOR AND DATASET 

A detailed description of the BABAR detector and of the 
algorithms used for charged and neutral particle recon- 
struction and identification is provided elsewhere [20, 21]. 
Charged particles are reconstructed by matching hits 
in the 5-layer double-sided silicon vertex tracker (SVT) 



IV. SIGNAL RECONSTRUCTION 

We reconstruct D° — > K~e + v e (^) decays in e + e~ — > 
cc events where the D° originates from the D* + — > D°w + . 
The main sources of background arise from events with 
a kaon and electron candidate. Such events come from 
T(4S I ) decays and the continuum production of charmed 
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hadrons. Their contribution is reduced using variables 
sensitive to the particle production characteristics that 
are different for signal and background events. 



A. Signal selection 



Charged and neutral particles are boosted to the cen- 
ter of mass system (cm.) and the event thrust axis is 
determined. The direction of this axis is required to be 
in the interval | cos(# t hrust)| < 0.6 to minimize the loss of 
particles in regions close to the beam axis. A plane per- 
pendicular to the thrust axis is used to define two hemi- 
spheres, equivalent to the two jets produced by quark 
fragmentation. In each hemisphere, we search for pairs 
of oppositely charged leptons and kaons. For the charged 
lcpton candidates we consider only electrons or positrons 
with cm. momentum greater than 0.5 GeV/c. 

Since the v e momentum is unmeasured, a kinematic 
fit is performed, constraining the invariant mass of the 
candidate e + K~v e system to the D° mass. In this fit, 
the D° momentum and the neutrino energy are estimated 
from the other particles measured in the event. The D° 
direction is taken as the direction opposite to the sum of 
the momenta of all reconstructed particles in the event, 
except for the kaon and the positron associated with the 
signal candidate. The energy of the jet is determined 
from the total cm. energy and from the measured masses 
of the two jets. The neutrino energy is estimated as the 
difference between the total energy of the jet and the 
sum of the energies of all reconstructed particles in the 
hemisphere. A correction, which depends on the value 
of the missing energy measured in the opposite jet, is 
applied to account for the presence of missing energy due 
to particles escaping detection, even in the absence of a 
neutrino from the D° decay. 

The D° candidate is retained if the x 2 probability of 
the kinematic fit exceeds 10 -3 . Detector performance for 
the reconstruction of the D° direction and for the missing 
energy are measured using events in which the D° decays 
into K~ir + . Corrections are applied to account for ob- 
served differences between data and simulation. Each 
D° candidate is combined with a charged pion, with the 
same charge as the lepton, and situated in the same hemi- 
sphere. The mass difference S(m) = m(D°ir + ) — m(D°) 
is evaluated and is shown in Fig. 1 . This distribution con- 
tains events which in addition pass the requirements on 
the Fisher discriminant F BB suppressing BB background 
and also give a satisfactory kinematic fit constraining the 
invariant mass . This last requirement is the reason of 
the slow decrease of the S(m) distribution. At large 5(m) 
values, a small excess of background is measured in data 
and the simulation is rescaled accordingly. Only events 
with 6(m) < 0.16 GeV/c 2 are used in the analysis. 
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FIG. 1: Comparison of the S(m) distributions from data and 
simulated events. MC events have been normalized to the 
sample luminosity according to the different cross sections. 
An excess of background events of the order of 5% is observed 
for large values of S(m). The arrow indicates the additional 
selection applied for the q 2 distribution measurement. 



B. Background rejection 

Background events arise from T(4S) decays and 
hadronic events from the continuum. Three variables 
are used to reduce the contribution from BB events: 
R,2 (the ratio between the second and zeroth order Fox- 
Wolfram moments [26]), the total charged and neutral 
multiplicity and the momentum of the soft pion (tt s ) from 
the D* + . These variables exploit the topological differ- 
ences between events with B decays and events with cc 
fragmentation. The particle distribution in T(45) decay 
events tends to be isotropic as the B mesons are pro- 
duced near threshold, while the distribution in cc events 
is jet-like as the cm. energy is well above the charm 
threshold. This also results in a softer D* + momentum 
spectrum in T(AS) decays compared to cc events. 

Corresponding distributions of these variables for sig- 
nal and background events are given in Fig. 2. These 
variables have been combined linearly in a Fisher dis- 
criminant. The requirement F BB > 0.5 retains 65% of 
signal and 6% of BB- background events. 
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FIG. 2: MC simulations of distributions of the variables used 
in the Fisher discriminant analysis to reduce the BB event 
background: a) the normalized second Fox- Wolfram moment 
(R2), b) the event particle multiplicity, c) the slow-pion mo- 
mentum distribution, in the cm. frame, d) the Fisher variable 
for BB and for charm signal events. 

Background events from the continuum arise mainly 
from charm particles as requiring an electron and a kaon 
reduces the contribution from light-quark flavors to a low 
level. Because charm hadrons take a large fraction of the 
charm quark energy charm decay products have higher 
average energies and different angular distributions (rel- 
ative to the thrust axis or to the D direction) compared 
with other particles in the hemisphere emitted from the 
hadronization of the c and c quarks. These other parti- 
cles are referred to as "spectator" in the following; the 
"leading" particle is the one with the largest momentum. 
To reduce background from cc events, the following vari- 
ables are used: 

• the D° momentum; 

• the spectator system mass, m sp ., which has lower 
values for signal events; 

• the direction of the spectator system momentum 
relative to the thrust axis cos # sp ._ thrust! 

• the momentum of the leading spectator track; 

• the direction of the leading spectator track relative 
to the D° direction; 



FIG. 3: MC simulations of some of the variables used in 
the Fisher discriminant analysis to reduce the cc-event back- 
ground: a) the D° momentum after the kinematic fit, b) the 
mass of the spectator system (peaks, at low mass values cor- 
respond to events with a single charged pion or photon recon- 
structed in the spectator system), c) the cosine of the angle 
between the spectator system momentum and the thrust di- 
rection, d) the cosine of the angle of the positron direction, 
relative to the kaon direction, in the ev e cm. frame. 

• the direction of the leading spectator track relative 
to the thrust axis; 

• the direction of the lepton relative to the kaon di- 
rection, in the dilepton rest frame, cos 9 e ; 

• the charged lepton momentum, p ei in the cm. 
frame. 

The first six variables depend on the properties of c-quark 
hadronization whereas the last two are related to de- 
cay characteristics of the signal. Distributions for four 
of the most discriminating variables are given in Fig. 3. 
D° — > K~ir + events have been used to tune the simu- 
lation parameters so that distributions of the variables 
used to reject background agree with those measured 
with data events. These eight variables have been com- 
bined linearly into a Fisher discriminant variable (F cc ) 
and events have been kept for values above 0. This selec- 
tion retains 77% of signal events that were kept by the 
previous selection requirement and rejecting 66% of the 
background (Fig. 4). 

The remaining background from cc-events can be di- 
vided into peaking (60%) and non-peaking (40%) candi- 
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FIG. 4: Distribution of the values of the Fisher variable in 
the signal region (8(m) < 0.16 GeV/c 2 in a), and for masses 
above the signal region (5(m) > 0.16 GeV/c 2 in b). 




q 2 r (GeV 2 ) 

FIG. 5: The measured g 2 . distribution (data points) compared 
to the sum of the estimated background and of the fitted 
signal components. 



dates. Peaking events are those background events whose 
distribution is peaked around the signal region. These 
are mainly events with a real D* + in which the slow 7r + 
is included in the candidate track combination. Back- 
grounds from e + e~ annihilations into light uu, dd, ss 
quarks and BB events are non-peaking. These compo- 
nents, from simulation, are displayed in Fig. 1. 



C. q 2 measurement 

To improve the accuracy of the reconstructed D° mo- 
mentum, the nominal D* + mass is added as a constraint 
in the previous fit and only events with a \ 2 probabil- 
ity higher than 1% are kept (Fig. 1 is obtained requiring 
only that the fit has converged). The measured q 2 distri- 
bution, where q 2 — [pr> — Pk) 2 , is given in Fig. 5. There 
are 85260 selected D° candidates containing an estimated 
number of 11280 background events. The non-peaking 
component comprises 54% of the background. 

To obtain the true q 2 distribution, the measured one 
has to be corrected for selection efficiency and detector 
resolution effects. This is done using an unfolding algo- 
rithm based on MC simulation of these effects. 

The variation of the selection efficiency as a func- 
tion of q 2 is given in Fig. 6. The resolution of the 
q 2 measurement for signal events is obtained from MC 



simulation. The resolution function can fitted by the 
sum of two Gaussian functions, with standard deviations 
(7i = 0.066 GeV 2 and cr 2 = 0.219 GeV 2 , respectively. 
The narrow component corresponds to 40% of the events. 

To obtain the unfolded q 2 distribution for signal events, 
corrected for resolution and acceptance effects, the Singu- 
lar Value Decomposition (SVD) [27] of the resolution ma- 
trix has been used. This method uses a two-dimensional 
matrix which relates the generated q 2 distribution to the 
detected distribution, q 2 , as input. After subtracting the 
estimated background contribution, the measured binned 
q 2 distribution is linearly transformed into the unfolded 
q 2 distribution. This approach provides the full covari- 
ance matrix for the bin contents of the unfolded distri- 
bution. Singular values (SV) are ordered by decreasing 
values. These values contain the information needed to 
transform the measured distribution into the unfolded 
spectrum, along with statistical uncertainties from fluc- 
tuations. Not all SV are relevant; non-significant values 
have zero mean and standard deviation equal to unity 
[27]. Using toy simulations, we find that seven SV have 
to be kept with events distributed over ten bins. Because 
the measurement of the form-factor parameters relies on 
the measured q 2 distribution, it does not require unfold- 
ing, and is independent of this particular choice. 
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FIG. 6: The efficiency as a function of q 2 , measured with 
simulated signal events, after all selection criteria applied. 



V. RESULTS ON THE Q 2 DEPENDENCE OF 
THE HADRONIC FORM FACTOR 



FIG. 7: Comparison between the normalized unfolded q dis- 
tribution obtained from this analysis, and those corresponding 
to two fitted models. Lower plot gives the difference between 
measured and fitted distributions. The error bars represent 
statistical errors only. 



The unfolded q distribution, normalized to unity, is 
presented in Fig. 7 and in Table II. Also given in this 
table are the statistical and total uncertainties and the 
correlations of the data in the ten bins. Figure 7 shows 
the result of fits to the data for two parameterizations of 
the form factor with a single free parameter, the simple 
pole and the modified pole ansatz. Both fitted distribu- 
tions agree well with the data. 

A summary of these and other form factor parame- 
terizations is given in Table III. These results will be 
discussed in detail in Section VB. 

The fit to a model is done by comparing the number 
of events measured in a given bin of q 2 with the expecta- 
tion from the exact analytic integration of the expression 
|p"R-(g 2 )| 3 |/ + (g 2 )| 2 over the bin range, with the overall 
normalization left free. The result of the fit correspond- 
ing to the parameterization of the form factor using two 
parameters (see Eq. 10) is given in Fig. 8. 




pole 



A. Systematic Uncertainties 

Systematic uncertainties of the form factor parame- 
ters are likely to originate from imperfect simulation of 
c-quark fragmentation and the detector response, from 
uncertainties in the background composition and the in- 
dividual contributions for the selected signal sample, the 



FIG. 8: Contours at 70% and 90% CL resulting from the fit of 
the parameterization of the form factor q 2 dependence with 
two parameters as given in Eq. (10). The value <5 po i c = 
corresponds to the modified pole ansatz. 



uncertainty in the modeling of the signal decay and the 
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TABLE II: Statistical and total uncertainty matrices for the normalized decay distribution (corrected for acceptance and finite 
resolution effects) in ten bins of q 2 from to 2 GeV 2 , and for the ratio Rd (see Section VI). The total decay distribution 
has been normalized to unity for q 2 varying over ten 0.2 GeV 2 intervals. The uncertainty matrices are provided for both the 
statistical (upper half) and total (lower half) uncertainties. The uncertainty on each measured value is given along the diagonal. 
Off-diagonal terms correspond to the correlation coefficients. 
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TABLE III: Fitted values of the parameters corresponding to different parameterizations of f+(q ). The last column gives the 
X 2 /NDF of the fit when using the value expected for the parameter. 



Theoretical Unit Parameters x^/NDF Expectations 

ansatz [x 2 /NDF] 

z expansion n = -2.5 ± 0.2 ± 0.2 5.9/7 

r 2 = 0.6 ±6. ±5. 
Modified pole a polc = 0.377 ± 0.023 ± 0.029 6.0/8 

Simple pole GeV/c 2 m po i c = 1.884 ± 0.012 ± 0.015 7.4/8 2.112(243/9] 
ISGW2 GeW 2 a/ = 0.226 ±0.005 ±0.006 6.4/8 0.104 [800/9] 



measurement of the q 2 distribution. We study the origin 
and size of various systematic effects, correct the MC sim- 
ulation, if possible, assess the impact of the uncertainty 
of the size of correction on the fit results, and adopt the 
observed change as a contribution to the systematic un- 
certainty on the fitted parameters for the different param- 
eterizations under study. Some of these studies make use 
of standard BABAR evaluations of detection efficiencies, 
others rely on special data control samples, for instance 
hadronic decays D° — ► K~n + or K^tt + -k°. 



1. c-quark hadronization tuning 

The signal selection is based on variables related to 
c-quark fragmentation and decay properties of signal 
events. Simulated events have been weighted to agree 
with the distributions observed in data. Weights have 
been obtained using events with a reconstructed D° de- 
caying into K~ir + . After applying these corrections, 
the distribution of the Fisher discriminant that contains 
these variables is compared for data and simulation. The 
remaining difference is used to evaluate the correspond- 
ing systematic uncertainty, ft corresponds to the varia- 
tions on fitted quantities obtained by correcting or not 
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for this difference, which is below 5% over the range of 
this variable. 



2. Reconstruction algorithm 

It is important to verify that the q 2 variation of the 
selection efficiency is well described by the simulation. 
This is done by analyzing D° — » K^tt + t: as if they were 
K~e + v e events. The two photons from the ir° are re- 
moved and events are reconstructed using the algorithm 
applied to the semileptonic D° decay. The "missing" n° 
and the charged pion play, respectively, the roles of the 
neutrino and the electron. To preserve the correct kine- 
matic limits, it is necessary to take into account that 
the "fake" neutrino has the ir° mass and that the "fake" 
electron has the ir + mass. 

Data and simulated events, which satisfy the same 
analysis selection criteria as for Kev e , have been com- 
pared. For this test, the cos(# e ) and p e are removed from 
the Fisher discriminant, because distributions for these 
two variables are different from the signal events. 

The ratio of efficiencies measured in data and simula- 
tion is fit with a linear expression in q 2 . The correspond- 
ing slope, (0.71 ±0.68)%, indicates that there is no signif- 
icant bias when the event selection criteria are applied. 
The measured slope is used to define a correction and to 
estimate the corresponding systematic uncertainty. 



3. Resolution on q 2 

To measure possible differences between data and 
simulation on the q 2 reconstruction accuracy, D° — > 
K ~7T + ir° events are used again. Distributions of the dif- 
ference q 2 — q 2 , obtained by selecting events in a given 
bin of q 2 are compared. These distributions are system- 
atically slightly narrower for simulated events and the 
fraction of events in the distant tails are higher for data 
(see Fig. 9). 

With the D° — » K ~ir + sample we study, in data and 
simulation, the accuracy of the D° direction and missing 
energy reconstruction for the D° — > K~e + v e analysis. 
This information is used in the mass-constrained fits and 
thus influences the q 2 reconstruction. Once the simula- 
tion is tuned to reproduce the results obtained on data 
for these parameters, the q 2 resolution distributions agree 
very well, as shown in Fig. 9. One half of the measured 
variation on the fitted parameters from these corrections 
has been taken as a systematic uncertainty. 



4- Particle identification 

Effects from a momentum-dependent difference be- 
tween data and simulated events on the charged lepton 
and on the kaon identification have been evaluated. Such 
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FIG. 9: Distribution of the difference between the true and 
the reconstructed value of q 2 . D° — > K~tt + tv° data events 
correspond to dark squares and open circles are used for 
simulated ones. Ratio (Data/MC) of the two distributions 
are displayed. The distributions in a) compare data and 
simulated events before applying corrections measured with 
D° — * K~tv + events, whereas these corrections have been ap- 
plied for plots in b). 
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differences, which are typically below 2%, have been mea- 
sured for selected, high purity samples of electrons and 
kaons. These corrections have been applied and the ob- 
served variation has been taken as the estimate of the 
systematic uncertainty. 

5. Background estimate 

The background under the D* + signal has two compo- 
nents that have, respectively, non-peaking and peaking 
behavior. 

The non-peaking background originates from non-cc 
events and from continuum charm events in which the 
tt s candidate does not come from a cascading D* + . By 
comparing data and simulated event rates for 5{m) > 
0.18 GeV/c 2 (see Fig. 1), a correction of 1.05 is deter- 
mined from simulation for the non-peaking background. 
This correction is applied and an uncertainty of ±0.05 is 
used as the corresponding systematic uncertainty. 

Events which include a slow pion originating from D* + 
decay contribute in several ways to the peaking back- 
ground. The production rate of D* + mesons in the sim- 
ulation is in agreement with expectations based on mea- 
surements from CLEO [28]. The uncertainty of ±0.06 on 
this comparison is dominated by the systematic uncer- 
tainty from the CLEO result. 

To study the remaining effects, the peaking back- 
ground components have been divided according to the 
process from which they originate and have been ordered 
by decreasing level of importance: 

• the K~ and the electron originate from a D° de- 
cay (54%). The main source comes from D° — > 
K~ir°e + v e . We have corrected the decay branch- 
ing fraction used for this channel in the MC (2.02%) 
using recent measurements (2.17±0.16% [1]). The 
uncertainty on this value has been used to evaluate 
the corresponding systematic uncertainty. 

• the electron comes from a converted photon or a 
Dalitz decay (24%). It has been assumed that the 
simulation correctly accounts for this component; 

• the K~ does not originate from a charm hadron 
(14%> ). This happens usually when there is another 
negative charged kaon accompanying the D* + . We 
have studied the production of charged kaon accom- 
panying a D* + using D° — > K~ir + events and mea- 
sure a correction factor of 0.87±0.02 and 0.53±0.02, 
respectively, for same sign and opposite sign K- 
D* pairs. The simulation is modified accordingly 
and the remaining systematic uncertainty from this 
source becomes negligible; 

• fake kaon candidate (mainly pions) (6%) or fake 
electrons (1%>). Differences between data and MC 
on the evaluation of fake rates have been studied 
in BABAR. As this affects small components of the 



total peaking background rate, the effect of these 
differences has been neglected. 



6. Fitting procedure and radiative events 

To fit form factor parameters we compare the num- 
ber of expected events in each bin with the measured 
one after all corrections. In this approach it is always 
assumed that the q 2 variation of f+(q 2 ) is given ex- 
actly by the form factor parameterization. This hy- 
pothesis is not correct, a priori, for radiative decays as 
q 2 = (pd-Pk) 2 = (Pe+Pv +P- ( ) 2 is not (perhaps) equal 
to the variable that enters in / + for such decays. PHO- 
TOS is used to generate decays with additional photons 
and the modified pole ansatz is taken to parameterize 
the hadronic form factor in signal events. To quantify 
possible distortion of the fit we compare the fitted value 
of a form factor parameter with the one obtained from a 
fit to the generated q 2 distribution (see Table IV). 

TABLE IV: Measured differences between the nominal and 
fitted values of parameters. Quoted uncertainties correspond 
to MC statistics. The last column gives the impact of the 
radiative effects on the form factor measurements as predicted 
by PHOTOS. 



Parameter 


Measured difference 


Bias from 




(true-fitted) 


radiation 


5(n) (xO.01) 


1.2 ± 13.0 


-0.4 + 2.7 


5(r 2 ) 


+ 1.7 ±3.6 


+ 1.9 + 0.7 


5(a po ie) (xO.01) 
S(m polc ) (McV/c 2 ) 
6{ai) (xO.OOlGeV" 2 ) 


-1.2 ± 1.4 


-1.1 + 0.3 


+4.5 ±6.3 
-2.7 + 3.1 


+4.6 + 1.4 
-2.4 + 0.7 



Corresponding corrections, given in the second column 
of Table IV, have been applied and quoted uncertainties 
enter in the systematic uncertainty evaluation. 

To evaluate the importance of corrections induced by 
radiative effects, we have compared also the fitted value 
of a parameter on q 2 distributions generated with and 
without using PHOTOS. Measured differences arc given 
in the last column of Tabic IV. They have not been 
applied to the values quoted in Table III for the different 
parameters. We measure also that radiative effects affect 
mainly the fraction of the decay spectrum in the first bin 
in Table II which has to be increased by 0.0012 to correct 
for this effect. 



7. Control of the statistical accuracy in the SVD approach 

Once the number of SV is fixed, one must verify that 
the statistical precision obtained for each binned un- 
folded value is correct and if biases generated by remov- 
ing information are under control. These studies are done 
with toy simulations. One observes that the uncertainty 
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TABLE V: Summary of systematic uncertainties on the fitted parameters. 



Source S (m po ie) <5a po i e $ («/) <5 (n) 
(MeV/c 2 ) (xO.Of) (xO.OOf GeV~ 2 ) (xO.Of) 


<5(r 2 ) 


c-hadronization tuning 


3.0 


0.6 


f.f 


6.7 


1.4 


Reconstruction algorithm 


7.8 


f.6 


3.f 


6.5 


0.2 


Resolution on q 2 


3.4 


0.7 


f.4 


3.0 


2.1 


Particle ID 


5.5 


IT 


2.3 


3.8 


0.4 


Background estimate 


7.7 


1.4 


3.0 


12.7 


2.0 


Fitting procedure 


6.3 


1.4 


3.1 


13.0 


3.6 


Total 


f5.0 


2.9 


6.0 


21.0 


4.9 



obtained from a fit of the unfolded distribution is un- 
derestimated by a factor which depends on the statistics 
of simulated events and is <~ 1.06 in the present anal- 
ysis. Pull distributions indicate also that the unfolded 
values, in each bin, have biases which are below 10% of 
the statistical uncertainty. Similar studies are done for 
the determination of form factor parameters. 



8. Summary of systematic errors 

The systematic uncertainties for determining form fac- 
tor parameters are summarized in Table V. 

The systematic error matrix for the ten unfolded val- 
ues is computed by considering, in turn, each source of 
uncertainty and by measuring the variation, Si, of the cor- 
responding unfolded value in each bin (i). The elements 
of the uncertainty matrix are the sum, over all sources 
of systematic uncertainty, of the quantities Si ■ Sj . The 
total error matrix is evaluated as the sum of the matrices 
corresponding, respectively, to statistical and systematic 
uncertainties. 



pole at higher mass (see Eq. 10) the two parameters are 
fully correlated and there is no unique solution, as illus- 
trated in Fig. 8. The modified pole ansatz corresponds 
to <5 po io = 0. 

In Table VI the fitted parameters for the simple pole 
ansatz and the modified pole [11] ansatz are compared 
for different experiments. The fitted pole masses are all 
well below the mass of the D* meson. The results pre- 
sented here are consistent within the stated uncertainties 
with earlier measurements. Except for the BELLE mea- 
surement, all other measurements appear to favor a value 
of Opoie that is lower than the value predicted by lattice 
QCD, namely a polc = 0.50 ± 0.04. 

TABLE VI: Fitted values for the parameters corresponding 
respectively to a pole mass and a modified pole mass model 
for the form factor. 



Experiment m po i c (GeV/c 2 ) Qpoic 
CLEO [4] 1.89 ± 0.05+o ot 0.36 ± 0.10+o '<£ 

FOCUS [5] 1.93 ± 0.05 ± 0.03 0.28 ± 0.08 ± 0.07 
BELLE [6] 1.82 ±0.04 ±0.03 0.52 ± 0.08 ± 0.06 
This analysis 1.884 ± 0.012 ± 0.015 0.38 ± 0.02 ± 0.03 



B. Comparison with expectations and with other 
measurements 

The summary of the fits to the normalized q 2 distri- 
butions are presented in Table III. As long as we allow 
the form factor parameters to be free in the fit, the fit- 
ted distributions agree well with the data and it is not 
possible to reject any of the parameterizations. 

However, if the form factor parameters are constrained 
to specific predicted values, the agreement is not good. 
For the ISGW2 model, the predicted dependence of the 
form factor on q 2 disagrees with the data and the fit- 
ted value of the parameter ai differs from the predicted 
value, ai = 0.104 GeV~ 2 by more than a factor two. 

As observed by previous experiments, the simple pole 
model ansatz, with m po i = mr)* = 2.112 GeV/c 2 does 
not reproduce the measurements. This means that the 
contribution from the continuum DK interaction can- 
not be neglected. If one introduces a second parame- 
ter (Jpoio to account for contributions from an effective 



In Fig. 10, the dependence of the form factor on q 2 is 
presented. The data are compared to earlier measure- 
ments by the FOCUS experiment, as well as with pre- 
dictions from lattice QCD calculations [19]. As stated 
above, the data favor a somewhat lower value for a po i e - 

The data have also been mapped into the variable z. 
Figure 1 1 shows the product P x $ x / + as a function of 
z. By convention, this quantity is constrained to unity 
at z = z max , which corresponds to q 2 = 0. We perform a 
fit to a polynomial, Px$x/ + ~Hnz + r 2 z 2 . The 
data are compatible with a linear dependence, which is 
fully consistent with the modified pole ansatz for f + (q 2 ), 
as illustrated in Fig. 11. 

VI. BRANCHING FRACTION MEASUREMENT 

The D — > K~e + v e branching fraction is measured 
relative to the reference decay channel, D° — > K ~n + . 
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FIG. 10: Comparison of the measured variation of 
/+(g 2 )//+(0) obtained in the present analysis and in the FO- 
CUS experiment [5]. The band corresponds to lattice QCD 
[19] with the estimated uncertainty. 



FIG. 11: Measured values for Px$x/+ are plotted versus —z 
and requiring that P x $ x /+ = 1 for z — z max . The straight 
lines represent the result for the modified pole ansatz, the fit 
in the center and the statistical and total uncertainty. 



Specifically, we compare the ratio of rates for the decay 
chains D*+ -> D°ir+,D° -> R-e + v e , and D° -> K~ir+ 
in data and simulated events, this way, many systematic 
uncertainties cancel, 



The third line corresponds to the ratios of measured 
numbers of signal events in data and in simulation, and 
the last line gives the ratios of the efficiencies to data and 
simulation. 



_ BR(D° 
D " BR(D Q 



■ K-e+u e ) data = BR(D° -> K- e + u e ) MC 
■* K-n+) d ^ BR(D° -> K-ty+)mc 
N(cc)Kev C(data)Kn 



,0 



N(cc)ktt C{data)Ke.v 

K-e+ Ve ) d ^ ^ N(D° -> A-7r+) MC 
N(D" -> /v-7r+) data 



x 



N(D° K~e+u e ) M c 

e(D° ->K-e + v e ) M c <D 
-T~F7\ _l \ x — 7~F\ 



(15) 



The first line in this expression is the ratio of the 
branching fraction for the two channels used in the sim- 
ulation: 



BR{D° -> K-e + v e ) MC _ 0.0364 
BR(D° -> K-tt+)mg ~ 0.0383' 



(16) 



The second line is the ratio of the number of cc simu- 
lated events and the integrated luminosities for the two 
channels: 



N{cc) Kev £(data) Kw 117.0 x 10 6 73.431b" 1 



N{cc) Kn C{dsXa) Kev 117.3 x 10 6 74.27 fb" 



r(i7) 



A. Selection of candidate signal events 

The selection of 13° — * K~e + v e candidates is explained 
in section IV A. For the rate measurement, the con- 
straint on the D* + mass is not applied and also the 
momentum of the soft pion candidate is not included 
in the Fisher discriminant variable designed to suppress 
SB-background. Since generic simulated signal events 
used in this measurement have been generated with the 
ISGW2 model, they have been weighted so that their q 2 
distribution agrees with the measurement presented in 
this paper. Furthermore, we require for the Fisher dis- 
criminant F B g > and restrict S(m) < 0.16 GeV/c 2 . 
After background subtraction, there remain 76283 ± 323 
and 95302 ± 309 events in data and simulation, respec- 
tively. This gives: 



N(D° -> K-e+v e ) data 
N(D° -> K-e+v e ) MC 



= 0.8004 ±0.0043. 



(18) 



To select D° — > K~n + candidates, the same data sam- 
ples are used and particles, in each event, are selected in 
the same way. The same selection criteria on the Fisher 
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FIG. 12: Events selected for the reference channel D* + — > D°ty + , D° — » .ft' 7r + . a) A'tt mass distribution for events selected in 
the range 5(m) E [0.143, 0.148] GeV/c 2 . b) S(m) distribution for events selected in the range m(KTv) £ [1.83, 1.89] GeV/c 2 . 
c) same distribution as in b), displayed on a larger mass range with the non-peaking background indicated (shaded area), d) 
S(m) distribution after non-peaking background subtraction, data (points with statistical errors) and simulated events (shaded 
histogram) . 



discriminant to suppress Bi3-events, on the thrust axis 
direction and on other common variables are applied. 
Events are also analyzed in the same way, with two hemi- 
spheres defined from the thrust axis. In each hemisphere 
a D° candidate is reconstructed by combining a charged 
K with a pion of opposite sign. These tracks have to 
form a vertex and the Kn mass must be within the range 
[1.77, 1.95] GeV/c 2 . Another charged pion of appropri- 
ate sign is added to form a D* + candidate. 

In addition, the following selection criteria are used: 

• the fraction of the beam momentum, in the cm. 
frame, taken by the D* candidate must exceed 0.48 
to remove contributions from BB events; 



• the measured D° mass must be in the range be- 
tween 1.83 and 1.89 GeV/c 2 . This requirement 
eliminates possible contributions from remaining 
D° — > K~K + or 7r + 7r~ decays (see Fig. 12-a); 

• the vertex fit for the D° and D* have to converge. 

The <5 (to) distribution for candidate events is shown in 
Fig. 12-c. The following components contribute to the 
D*+ signal (see Fig. 12-b): 

• D° — > K~ir + with no extra photon; 

• D° — > K~ir + with at least one extra photon; 

• D — > K~ir + (j) where the ir + , mainly from the 
D* + , decays into a muon. 
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The 6(m) distribution corresponding to other event 
categories does not show a peaking component in the 
D* + signal region. The total background level, is nor- 
malized to data using events in the S(m) interval between 
0.165 and 0.200 GeV/c 2 (see Fig. 12-c). This global scal- 
ing factor is equal to 1.069±0.011. After background sub- 
traction, the 5(m) distributions obtained in data and sim- 
ulation can be compared in Fig. 12-d. Since the D* + sig- 
nal is narrower in the simulation, we use a mass window 
such that this difference has a small effect on the mea- 
sured number of D* + events in data and in the simula- 
tion. There are 166960±409 and 134537±374 candidates 
selected in the interval 5{m) <E [0.142, 0.149] GeV/c 2 for 
simulated and data events respectively. This means: 

N(D° -> K-tt+)mc 



N(D° -> K-TT+) 



= 1.230 ±0.0046. 



data 



(19) 



B. Efficiency corrections 

The impact of the selection requirement on the recon- 
structed Kit mass, has been studied. The Kir mass dis- 
tribution signal for simulated events (not including radia- 
tive photons) is compared with the corresponding distri- 
bution obtained with data after background subtraction. 
The background contributions are taken from the simula- 
tion. The fraction of D° candidates in the selected mass 
range (between 1.83 and 1.89 GeV/c 2 ) is (97.64±0.25)% 
in MC and (97.13 ± 0.29)% in data events. The ratio of 
efficiencies is equal to: 

e(D° -> K-Ti+) data 



K- 



= 0.9947 ±0.0039. 



MC 
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e(D° 

Since D° — > K~e + v e events have been selected using 
a selection requirement on <5(m), we need to confirm that 
the distribution of this variable is similar in data and 
simulation. This is checked by comparing the distribu- 
tions obtained with D° — > K~ir + ir° events analyzed as 
if they were semileptonic decays. The 5(m) distributions 
are compared in Fig. 13. Below 0.16 GeV/c 2 , there are 
0.93552±0.00066of the D*+ candidates in the simulation 
and 0.93219 ± 0.00078 for data. The corresponding ratio 
of efficiencies (MC/data) is equal to 1.0036 ± 0.0010. 

Using D° — ► K~ir + ir° events, we also measure the dif- 
ference between the fraction of events retained after the 
mass-constrained fits. Namely, it is 0.98038 ± 0.00037 in 
the simulation compared to 0.97438 ± 0.00049 in data. 
The relative efficiency (MC/data) for this selection is 
1.0062 ± 0.0006. Based on these two measured correc- 
tions the ratio of efficiencies are: 

e(D° -» K- 



1.0098 ±0.0011. 



(21) 



e(£>° -> K-e+ Ve ) data 

The quoted uncertainties, in this section, are of statistical 
origin and will be included in the statistical uncertainty 
on Rr>. Other differences between the two analyzed chan- 
nels are considered in the following section and contribute 
to systematic uncertainties. 
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FIG. 13: S(m) distribution for D° -> K~ 'n + ir events ana- 
lyzed as if they were semileptonic decays. Distributions have 
been normalized to unity; note that the bin size is not uni- 
form. The bottom plot shows the ratio of the two distribu- 
tions above. 



C. Systematic uncertainties on Rd 

A summary of the systematic uncertainties on Rjj are 
given in Table VII. They originate from selection criteria 
which are different for the two channels. Some of these 
uncertainties are the same as those already considered 
for the determination of the q 2 variation of /+ . 



1. Correlated systematic uncertainties 

Systematic uncertainties on the decay rate coming 
from effects that contribute in the measurement of the 
q 2 dependence of f+(q 2 ) are evaluated in Section VA 
and the full covariance matrix for the measurements of 
the number of D° — > K~e + v e signal events and the frac- 
tion of the decay spectrum fitted in each of the ten bins 
is determined. Among the sources of systematic uncer- 
tainties, listed in Table V, those corresponding to: 

• the reconstruction algorithm, 

• the tuning of the resolution on q 2 , 

• the corrections applied on electron identification, 
and 
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• the background normalization 



D. Decay rate measurement 



are taken as common sources. Corresponding relative 
uncertainties on Rrj are given in Table VII. 

Other systematic uncertainties contributing to the 
form factor measurement also affect the reference chan- 
nel and so their effects on R D cancel. They are related to 
the c-hadronization tuning and to the corrections applied 
on the kaon identification. 



2. Selection requirement on the Fisher discriminant 

The stability of the fraction of D° — > K~e + u e events 
selected in data and in simulation as a function of the 
Fisher discriminant, F cc , designed to suppress cc back- 
ground has been examined. This is done by comparing 
the distributions of this variable measured in data and in 
simulation as given in Fig. 4 for two selected intervals in 
5(m). 

The value corresponding to F cc > and for events se- 
lected in the range S(m) < 0.16 GeV/c 2 is used as the 
central result and half the difference between the mea- 
surements corresponding to F cc greater than —0.25 and ± 
0.25 is taken as systematic uncertainty. This range cor- 
responds to a relative change of 40% of the efficiency for 
signal events, and gives an uncertainty of ±0.0061 on the 
ratio of data and simulated signal candidates. 



3. D* + counting in D° -> A~tt 4 



Combining all measured fractions in Eq. (15), the mea- 
sured relative decay rate is: 



R D = 0.9269 ± 0.0072 ± 0.0119. 



(22) 



Using the world 
BR(D° -> K-TT+) 
K 

the last quoted uncertainty corresponds to the accuracy 
on BR(D° -> K--K+). 



average for the branching fraction 
: (3.80 ±0.07)% [1], gives BR(D° -> 
■■ + v e {l)) = (3-522 ± 0.027 ± 0.045 ± 0.065)%, where 



VII. SUMMARY 

The decay rate distribution for the channel D a — > 
K~ e + v e ("f) has been measured in ten bins in Table II. 
Several theoretical expectations for the variation of this 
form factor with q 2 have been considered and values for 
the corresponding parameters have been obtained (see 
Table III). The q 2 variation of the form factor can be 
parameterized with a single parameter using different ex- 
pressions. The ISGW2 model with expected values for 
the parameters is excluded, as is the pole mass parame- 
terization with m po ic = mu» . 

The value of the decay branching fraction has been also 
measured independent of a model. 

Combining these measurements, the value of the 
hadronic form factor is obtained: 



D* + candidates are selected in the range S(m) G 
[0.142, 0.149] GeV/c 2 . From the simulation it is ex- 
pected that the fraction of signal events outside this 
interval is equal to 1.4%. Even though the D* + sig- 
nal is slightly narrower in the simulation, there is not 
a large discrepancy in the tails. The fraction of signal 
events measured in the sidebands S(m) E [0.140, 0.142]® 
[0.149, 0.150] GeV/c 2 is 0.4% and 0.5%, respectively, for 
simulation and data. An uncertainty of ±0.004, cor- 
responding to 30% uncertainty on the total fraction of 
events outside the selected 5{m) interval is assumed. 



TABLE VII: Summary of systematic uncertainties on the rel- 
ative decay rate measurement. 



Source 


Relative variation 


Reconstruction algorithm 


±0.42% 


Resolution on q 2 


~ 


Electron ID 


±0.56% 


Background subtraction 


±0.63% 


Cut on Fisher variable 


±0.76% 


D*+ counting (D° -> K~it + ) 


±0.40% 


Total 


±1.27% 



/+(0) = 



u. 



/24tt 3 BR 

G 2 F r DO r 



(23) 



where BR 
fraction, t d o 



is the measured D° - 
(410.1 ±1.5) x 10" 



■> K e + v e branching 
15 s [1] is the L>° life- 

To 



time and / = \p K (q 2 ) \" \f+(q 2 )/ f+(0) \ dq 2 

account for the variation of the form factor within one 
bin, and in particular to extrapolate the result at q 2 = 0, 
the pole mass and the modified pole ansatze have been 
used; the corresponding values obtained for /+(0) differ 
by 0.002. Taking the average between these two values 
and including their difference in the systematic uncer- 
tainty, this gives 



/+(0) = 0.727 ± 0.007 ± 0.005 ± 0.007, 



(24) 



where the last quoted uncertainty corresponds to the ac- 
curacy on BR(D° — ► K~n + ), t d q and |U CS |. It agrees 
with expectations and in particular with LQCD compu- 
tations [19]. Using the z expansion of Eq. (6), we find 
a = (2.98 ± 0.01 ± 0.03 ± 0.03) x 10~ 2 . 

The high accuracy of the present measurement will be 
a reference test for improved lattice determinations of the 
q 2 variation of /+ . 
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